Conductivity of Doped Two-Leg Ladders 
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Recently, conductivity measurements were performed on the hole-doped two-leg ladder material 
Sri4_ x Ca x Cu24C>4i. In this work, we calculate the conductivity for doped two- leg ladders using a 
model of hole-pairs forming a strongly correlated liquid - a single component Luttinger liquid - 
in the presence of disorder. Quantum interference effects are handled using renormalization group 
methods. We find that our model can account for the low energy features of the experimental 
results. However, at higher energies the experiments show deviations from the predictions of this 
model. Using the results of our calculations as well as results on the ground state of doped two-leg 
ladders, we suggest a scenario to account for the higher energy features of the experimental results. 



Over the last few years Jr there has been considerable in- 
terest in ladder materials.El In particular, the two-leg lad- 
der has received quite a bit of attention. Theoretically, 
the two-leg ladder is popular becase one can do controlled 
calculations, both analytically and numerically. Experi- 
mentally, two-leg ladders are popular because there are 
several known two-leg ladder materials, and these mate- 
rials can be fabricated in a controlled way. In particular, 
the two-leg ladder material Sri4_ x Ca x Cu«404i has re- 
ceived a considerable amount of attention^ Unlike some 
of the other ladder materials which must be fabricated 
under high pressure, this material can be fabricated un- 
der ambient pressure. Also, when Ca replaces some of 
the Sr, holes are doped into the ladders. However, these 
holes were already present in the material in the charge 
reservoir layers composed of Cu02 chains; doping with 
Ca simply moves them to the ladders. Therefore, this 
material is self-doped, and hence, extremely clean. Be- 
sides these practical reasons, the two-leg ladder is inter- 
esting because its properties are remarkably similar to 
the cuprate superconductors, namely when doped with 
holes, the system has a spin gap and the doped holes 
form pairs with a d 2 .2_ y 2-like internal structure. 

Recently, conductivity measurements were per- 
formed on the hole-doped two-leg ladder material 
Sri4_ x Ca x Cu2404i Jj The c-axis resistivity, p c , (i.e. the 
resistivity parallel to the chains of the ladder) is insu- 
lating at low temperatures and shows metallic T depen- 
dence above To « 1QQK for x — 11; p c is nonmetallic 
over the entire region T < 300K for x = 8. The a-axis 
resistivity, p a , (i.e. the resistivity parallel to the rungs of 
the ladder) is insulating at low temperatures and shows a 
gradual change of sign in the T coefficient at T* ss 250X 
for x — 11; p a is insulating for x = 8. (See Fig. 1 of 
Ref. |.) 

The a-axis optical conductivity, cr± a (u>), shows pseudo- 
gap behavior in the low frequency region, the depressed 
spectral weight being transferred into the higher energy 
region. (See Fig. 2 of Ref. ||.) The c-axis optical con- 
ductivity, <7i c (w), shows a low frequency peak. Unlike a 
Drude peak, the peak in a\ c {uj) is located at a finite fre- 
quency, loq: loq ~ lOOcm -1 for x — 8 and ujq ~ 50cm -1 
for x = 11. (See Fig. 3 of Ref. |.) 



In this work, we calculate the conductivity for doped 
two-leg ladders using a model of hole-pairs forming a 
strongly correlated liquid. We focus on the c-axis con- 
ductivity, but we make some remarks about the a-axis 
conductivity. Our model is motivated by several known 
facts. First, the two-leg ladder (at half-filling and for 
not too large dopings) is known to have a spin gap. Sec- 
ond, the holes doped into the two-leg ladder are known to 
form pairs with an internal d x 2 _ y 2-\\k& structure. There- 
fore, for our (effective low-energy) model, we consider a 
liquid of d x 2_ y 2 hole pairs in the presence of (weak) dis- 
order. Due to the one-dimensional nature of the ladder, 
the pairs do not condense. Rather, they form a strongly 
correlated liquid — a single-component Luttinger liquid. 

It should be noted that our model is not purely phe- 
nomenological. At weak coupling, controlled calculations 
using bosonization have been doneou In these calcula- 
tions, one gets four bosonic modes — two spin modes 
and two charge modes. These calculations show that, 
for not too large dopings, both of the spin modes are 
gapped and one of the charge modes is gapped. There- 
fore, the low-energy physics is described byra single gap- 
less bosonic charge mode. Numerical workH shows that 
there is no phase transition as one goes from weak to 
strong coupling. Therefore, the low-energy effective the- 
ory in terms of a single gapless bosonic charge mode is 
still valid. However, the parameters of this effective the- 
ory cannot be computed as they can in weak coupling; 
they must be taken as phenomenological parameters from 
numerics. 

The Hamiltonian we consider is 
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This Hamiltonian was first considered in Ref. ^, and in 
that work they computed the resistivity, p c (T). In this 
work, we extend their calculation of p c (T) to lower tem- 
peratures, and we also calculate the optical conductivity, 
a\ c {uj). In Eq. 0, tj> and it are conjugate fields satisfying 
[4>{x), Tt(y)] = iSLx — y). £(x) and £*(a:) are Gaussian 
random variablesEI satisfying 
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£*(xi)£(x 2 ) = Ap S(xt - x 2 ) , 



(2) 



where the overbar denotes disorder averaging. K, the 
Luttinger liquid parameter, is determined by the inter- 
actions. For attractive interactions, K > 1; for repulsive 
interactions, K < 1. (-ft' = 1 for a non- interacting sys- 
tem.) Note that the coupling to disorder is via the 4fcp 
charge density wave (CDW).B Due to the gaps in the spin 
modes and one of the charge modes, the 2/cf CDW cou- 
pling is irrelevant; the most relevant coupling generated 
is to the 4/cf CDW. For later coniiepience, we define a 
dimensionless disorder strength, I?Ji3 



V = 



1 
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where rii is the density of impurities and Eq is a high- 
energy cut-off. Since our model is a liquid of d x 2_ y 2 hole 
pairs, Eq ~ pair binding energy. 

We begin with a standard renormalization group (RG) 
analysis. In an RG treatment we coarse-grain our sys- 
tem, integrating out high energy degrees of freedom, 
and deriving an effective theory for the low energy de- 
grees of freedom. The RG equations for the parameters 
can be obtained in a straightforward way using standard 
techniques .13 They are 
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dK 
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dv 
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A few words are in order about the RG equations. From 
the RG equation for V, we see that there is a localization- 
delocalization transition at K = 3/2. However, for the 
systems we are considering, K < 1. Hence, T> is strongly 
relevant and flows to strong coupling; all states will be 
localized at low temperatures. Also note that K as well 
as the speed of the excitations, v, decrease under the RG. 

From the Kubo formula, the c-axis conductivity for our 
model is given by 
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oi c (w, T) — — lm(J(x, t)J(x',t'))(q = 0, u) 
ijj 

where J(x, t) is the current operator 



J(x,t) = — = d t (, 

\/7T 



Fourier transforming in space and time, we have 



(5) 



<7i c (w,T) = Im(0(x,i)0(x',t'))(9 = 0,w) . (6) 

From the RG analysis, we saw that the disorder grows 
strong at low energies and the system localizes. In this 
work, we will not be concerned with the properties of 
the localized (insulating) state. More formally, we will 
be interested in the properties near the Luttinger liq- 
uid (free boson) fixed point. Therefore, we can use RG 



techniquesH£3 to compute the conductivity. Consider the 
Green's function 



G{w,V) = {<j>(x,t)4>(x',t')}(q = 0,u) 



Suppose we scale energies — ui 
Green's function transforms as 
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G(uj,V) = b 2 Z G{bcu,b 3 - 2K V) 



(7) 
The 
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The factor b 2 comes from dimensional analysis, while 
the factor Z is due to the "wave function renormaliza- 
tion" which arises in the RG. For simplicity, we ignore 
the renormalization of K and v. Since, by assump- 
tion, the disorder is weak and since we are far from 
the localization-delocalization transition, the renormal- 
ization of K and v should be small.El Writing Z = b 1 , 
where 7 = T>K, and iterating gives 



G(uj, V) = exp 



Jo 



dl (2-7(0) 
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We choose I* such that e l * = (E /E), where E 
T. 

Notice that on the right hand side of Eq. |^, G is evalu- 
ated at Eq. At such high energies, quantum interference 
effects should be negligible. In this regime, the Born ap- 
proximation gives accurate results for the conductivity. 
In the Born approximation, the Green's function has the 
form 



G(u,V) 



iuj/r 



(10) 



The functional form of r can be deduced from scaling 
arguments I 11 !" 2 ! In the Born approximation, 1/r ~ riiV T>. 
Also, from scaling 
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T(bE,b 3 ~ 2K V) t(E,V)' 
where E = lu or T. Therefore, we deduce that 
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Knowing the conductivity in the Born approximation, we 
can use Eq. ||to relate the conductivity at lower energies, 
where quantum interference effects are important, to the 
conductivity at higher energies, where quantum interfer- 
ence effects are negligible and the Born approximation is 
accurate. 

Let us first understand what our expressions tell us 
about the resistivity, p c {T). Using Eqs. 10 and [l2| with 
Eq. H gives 



Pc(T) 



2-2K 



exp 



KV 

; - 2K 



{l/xf~ 2K - 1 



(13) 
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where we have introduced x = T/Eq. Consider the first 
factor on the right hand side. ThiSr-aower-law factor can 
be deduced by scaling arguments. liil The second factor 
arises from the "wave function renormalization" in the 
RG and is due to quantum interference/localization ef- 
fects. 

Now let us consider the optical conductivity, aj c (u>) (at 
T = 0). Again, using Eqs. Eol and |l| with Eq. | gives 



<j 1c {uj) 
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exp 
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(l/x) 



3-2K 



- 1 



(14) 



where 8 ~ riiV T> / E and we have introduced x = w/E . 
Consider the first factor on the right hand side. No- 
tice that for K = 1, we have a peak at u> — — a 
Drude peak. For K < 1 (i.e. for an interacting system 
with repulsive interactions), the peak is shifted to finite 
frequency. The smaller the value of K, the higher the 
frequency at which the peak occurs. Now consider the 
second factor. This factor arises from the "wave function 
renormalization" in the RG. This factor is due to quan- 
tum interference/localization effects, and it suppresses 
the low frequency behavior. 
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FIG. 1. p(T) vs. T. Solid lines are from Eq. 
experimental data is from Ref. |^. 
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How do our results compare with the experimental re- 
sults of Ref. H? In Fig. we plot our calculated result 
for p c {T) as well as the experimental results of Ref. [|. 
To fit the data for x = 8, we took K = .25, V = 0.1, 
E = 425K, and p Q = 0.74mf2 cm. To fit the data for 
x = 11, we took K = 0.78, V = 0.1, E = U2K, and 
po = 0.285mSl cm. For x = 8, we were able to get a 
reasonable fit to the data at lower temperatures. How- 
ever, we cannot account for the rather constant resistivity 
starting at ~ 200K . We believe that pairing fluctuations 
and the renormalization of K and v, which we ignored 
in our calculation, could account for this difference. For 
x = 11, we were able to get a reasonable fit to the data at 
low temperatures. However, our calculated result cannot 



account for the metallic behavior seen above To « 100/G 
We believe the metallic behavior is due to pairing fluctu- 
ations, which were not included in our model. (See below 
for more details on pairing fluctuations.) 
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FIG. 2. ai(u) vs. Wave number. Solid lines are from 
Eq. [li]. The experimental data is from Ref. ^| at 10K . 



In Fig. g we plot our calculated result for ai c (uj) as 
well as the experimental results of Ref. ||. We see that 
<7i c (w) shows a narrow, finite frequency peak. This finite 
frequency (rather than u> = 0) peak is due to the strong 
correlations; the smaller the value of K, the larger the 
frequency at which the peak occurs. However, at higher 
energies, our calculated result does not reproduce some 
of the features of the measured conductivity. Namely, 
there is appreciable spectral weight at higher energies, 
while the calculated result gives practically no spectral 
weight at higher energies. We believe this appreciable 
spectral weight at higher energies is due to pairing fluc- 
tuations, which were not included in our model. In order 
to fit the data, we took K = 0.25, V = 0.1, 6 = 0.075, 
E a = 0.04cm -1 , and er = 150fl~ 1 cm^ 1 for x=8; we 
took K = 0.78, V = 0.1, 8 = 0.1, E = 0.035cm -1 , 
and Co — 355f2 -1 cm -1 for x=ll. The values for Eg and 
Co are inconsistent with the values determined from the 
resistivity. For x = 8, the peak occurs at ~ 100cm -1 
(= 144i-Q. At these energies the renormalization of K 
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and v, which we ignored in our calculation, could be com- 
ing into play. For x = 11, the peak occurs at ~ 50cm -1 
(= 72K). At these energies, pairing fluctuations are 
probably still be present and could account for our in- 
ability to fit the data consistently. (See below for more 
details on pairing fluctuations.) 

For not too large dopings, a simple picture which cap- 
tures the essential physics of the ground state of doped 
two-leg ladders is hole pairs moving through a nearest^ 
neighbor resonating valence bond (RVB) background .lL2l 
A typical configuration in the ground state is shown in 
Fig. In this picture, the spin gap is related to the en- 
ergy necessary to break the short ranged singlet bonds; 
pairing also occurs in order to minimize the number of 
broken singlet bonds. From this picture of the ground 
state, combined with our calculations for the conductiv- 
ity, we suggest a scenario to explain the experimental 
results of Ref. ||. 




FIG. 3. Typical configuration in the ground state of doped 
two-leg ladders - hole pairs mo ving through an RVB back- 
ground. (Figure taken from Ref. O .) 

To (« 100-ftT) is the temperature atmwhich the c-axis 
resistivity becomes metallic for x — 11.0 We believe that 
this is the temperature at which pairing fluctuations oc- 
cur. In other words, the RVB background is still, for 
the most part, in tact. However, the holes are no longer 
tightly bound into pairs. In the bosonization treatment, 
pairing is due to the "pinning" of p - = 9 p j, — pa , the 
relative phase of thfi-charge modes in the bonding and 
antibonding bandsmJ When 9 p - is pinned, the relative 
phase between the bonding and antibonding bands be- 
come locked together, with fluctuations being gapped. 
In our calculation, we took p - to be pinned; we ignored 
fluctuations, since they are gapped. We speculate that 
above To, 9 p - is no longer pinned and is free to fluctu- 
ate. This also explains why our calculated conductivity 
has almost no spectral weight at higher energies. Namely, 
at higher energies the missing spectral weight is due to 
the fluctuations of 9 p - which we ignored. 

T* (w 250K) is the temperature at which the a- 
axis resistivity becomes metallic for x = 11; it was 
found to coincide with the spin-gap inferred from NMR 
measurements.^ From the ground state wave function, wc 
see that it is energetically costly for single electrons to 
tunnel between neighboring ladders as long as the RVB 
background is in tact i.e. as long as there is a spin-gap. 
Above T* , the spin-gap is destroyed, and hence it is not 
costly for single electrons to tunnel between neighbor- 
ing ladders. In the bosonization language, due to the 
spin gap, single particle tunneling is irrelevant. How- 
ever, the operator which tunnels a pair between neigh- 
boring ladders is generated under the RG. It is the most 
relevant operator generated, and hence dominates the a- 
axis transport. Once the spin-gap is destroyed, the a-axis 



transport will be dominated by single particle tunneling. 

cia(w) was observed to exhibit pseudogap behavior. 
The pseudogap was identified to open up at ~ 1100cm -1 
(= 1584A) for x = 8 and - 600cm- 1 (= 864A) for 
x = 1 1 .0 These values are much larger than the spin-gap 
values. However, it is possible that the energies at .which 
the pseudogap begins were overestimated in Ref. ||.EJ We 
believe that the pseudogap should coincide with the spin- 
gap, namely because that is where single particles become 
free to tunnel between neighboring ladders. 

To summarize, we computed the conductivity of doped 
two-leg ladders. Our model consisted of hole pairs form- 
ing a single component Luttinger liquid in the presence 
of disorder; quantum interference effects were handled 
using RG methods. Our results could account for the 
low-energy features of Ref. || However, there were some 
discrepancies at higher energies, namely we could not ac- 
count for the metallic behavior in p c (T) for x = 11, and 
we could not account for the significant spectral weight in 
cti c (w) at higher energies. Based on the results of our cal- 
culations as well as results on the ground state of doped 
two-leg ladders, we suggested that these discrepancies 
were due to pairing fluctuations which were not included 
in our model. We also commented that it is quite nat- 
ural that T* corresponds to the spin-gap inferred from 
NMR measurements; single particles are free to tunnel 
between ladders above T*, but not below T*. Finally, we 
commented that we believe the energy at which the pseu- 
dogap begins in ai a (u>) should coincide with the spin-gap, 
and that its value may have been overestimated in Ref. |3J. 
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